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NATIONAL  ADVISORY  COMMITTEE  FOR  AERONAUTICS 


TECHNICAL  NOTE  NO.  1768 

A  RAPID  GRAPHICAL  METHOD  FOR  COMPUTING  THE  PRESSURE 
DISTRIBUTION  AT  SUPERSONIC  SPEEDS  ON  A  SLENDER 
ARBITRARY  BODY  OF  REVOLUTION 
By  Jim  Rogers  Thompson 

SUMMARY 


A  method  is  presented  by  which  the  pressure  distribution  on  a 
slender,  arbitrary  body  of  revolution  moving  at  supersonic  speeds  at 
zero  angle  of  attack  may  be  calculated  with  an  appreciable  saving  of 
time  as  compared  with  that  required  by  other  methods.  The  method  was  - 
developed  from  the  linearized  potential— flow  relations  given  by  Jones 
and  Margolis  in  NACA  TN  No.  108l  and  is  essentially  a  graphical  computing 
technique  based  on  the  assumption  that  a  body  may  be  represented  satis¬ 
factorily  by  a  relatively  small  number  of  point  sources  and  sinks .  An 
example  is  presented  to  illustrate  the  computational  procedure. 

The  method  gives  results  for  parabolic  bodies  in  excellent  agreement 
with  results  obtained  by  the  method  of  Jones  and  Margolis.  It  was  found 
that  a  parabolic  body  of  fineness  ratio  6  could  be  represented  satis¬ 
factorily  by  as  few  as  10  sources  and  sinks. 

Pressure  distributions  computed  for  the  Gorman  missile  A^VIP 
at  Mach  numbers  of  I.87  and  I.56  are  compared  with  wind-tunnel  measure¬ 
ments  and  with  results  of  calculations  by  the  characteristics  methods 
of  Sauer  and  Tollmein  to  provide  some  information  on  the  accuracy  and 
reliability  of  the  method  presented  herein  when  applied  to  arbitrary 
shapes.  Although  the  body  of  the  A4-V1P  is  somewhat  more  blunt  than 
would  be  desirable  for  treatment  by  linearized  methods,  the  method 
presented  reproduces  the  salient  features  of  the  measured  pressure 
distribution  with' comparatively  small  discrepancies. 


INTRODUCTION 


The  problem  of  computing  the  pressure  distribution  on  a  body  of 
revolution  moving  at  supersonic  speeds  has  been  treated  by  several 
authors  (references  1  to  6)  with  varying  degrees  of  approximation. 
References  1  to  5  present  methods  based  on  the  isentropic  linearized 
potential— flow  equation  first  applied  to  the  problem  by  Von  Karman 
in  reference  1.  A  more  rigorous  method  based  on  the  method  of  character¬ 
istics  is  also  presented  in  reference  5  and  is  generalized  by  Ferri 
in  reference  6  to  include  nonisentropic  rotational  flows.  The  usefulness 
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of  these  methods  for  the  design  of  supersonic  aircraft  or  missiles  has 
been  severely  limited,  however,  by  the  excessive  labor  involved  in 
the  computations. 

Jones  and  Margolis  (reference  4)  and  Lighthill  (reference  3)  have 
shown  that,  for  certain  mathematical  body  shapes,  the  pressure  distribution 
may  be  obtained  in  the  form  of  relatively  simple  equations.  In  general, 
however,  the  designer  cannot  utilize  these  simple  mathematical  shapes 
exclusively  because  of  other  requirements  of  the  design;  for  example, 
providing  vision  for  the  pilot  or  seeker,  tankage  space  equally  distributed 
about  the  center  of  gravity,  and  shape  changes  to  allow  simpler  or 
stronger  construction.  In  addition,  the  possibility  exists  that  body 
shapes  other  than  those  specified  by  simple  mathematical  relations  might 
have  the  lowest  drag,  particularly  when  these  shapes  are  combined  with 
wings  to  form  a  practical  airplane  or  missile  configuration.  The  need 
of  a  rapid  method  for  computing  the  pressure  distribution  on  arbitrary 
bodies  is  therefore  evident. 

An  approximate  method  for  computing  the  pressure  distribution  on 
slender,  arbitrary  bodies  of  revolution  moving  at  supersonic  speeds  and 
zero  angle  of  attack  Is  presented  herein.  Use  of  a  graphical  computing 
technique  allows  the  pressure  distribution  to  be  obtained  at  an  appreciable 
saving  in  time  as  compared  with  other  methods.  The  method  (which  is 
referred  to  hereinafter  as  the  "point-source  method")  was  developed  from 
the  results  of  Jones  and  Margolis  (reference  4)  and  is  based  on  a 
simplifying  assumption  concerning  the  source  distribution  used  to  represent 
the  body.  The  validity  of  this  assumption  in  the  case  of  parabolic  bodies 
is  examined  by  comparison  of  results  obtained  by  the  point-source  method 
with  results  obtained  by  the  method  of  Jones  and  Margolis. 

Although  results  obtained  by  the  point-source  method  are  subject  to 
the  uncertainties  inherent  in  the  linearized  potential— flow  equations 
and  to  additional  uncertainties  due  to  the  graphical  approximations  used 
in  the  computations,  the  method  is  believed  to  provide  a  quick  and  useful 
means  for  comparing  the  supersonic  characteristics  of  slender,  arbitrary 
bodies  of  revolution*)  Some  information  concerning  the  reliability  of 
the  results  was  obtained  by  comparison  of  pressure  distributions  computed 
by  the  point-source  method  with  calculations  by  the  characteristics  method 
and  with  wind-tunnel  measurements  presented  in  references  7  and  8  for  the 
German  missile  A4V1P. 

The  computational  procedure  for  the  point— source  method  is 
illustrated  by  an  example. 


SYMBOLS 

* 


P  pressure  coefficient 

x,r 


cylindrical  coordinates 


5 

fv(l) 

R 

F 

6 
Q 

Z^3C 

n 

k 

M 


abscissa  of  source  or  sink  (measured  along  x-axis) 

rate  of  change  of  source  Strength  (with  x) 

radius  of  body 

fineness  ratio  of  body 

half  angle  of  cone 

pressure  function  (f ’ (|)Ax) 

interred  length 

posit ire  integer 

constant 

Mach  number 


$  »  \/M2  -  1 

Subscripts  indicate  that  the  value  of  the  quantity  being  considered 
is  taken  at  that  particular  point .  Subscript  o  refers  to  the  nose  of 
the  body. 


ISpBLOlMiRff  OF  METHOD 


In  reference  4  Jones  and  Margblis  show  that  a  slender  body  of 
rerolution  moving  at  supersonic  speed  With  zero  angle  of  attack  may  be 
approximated  by  an  infinite  number  of  sources  or  sinks  located  on  the 
axis  of  the  body  and  that  the  pressure  coefficient  P  at  any  point  X 
on  the  body  surface  is  given  by 


where  ?*({)  is  the  first  deriratire  of  the  function 


Explicit  relations  may  be  easily  Obtained  from  equations  (l)  and  (2)  for 
the  pressure  distribution  on  cones  and  on  parabolic  bodies .  For  arbitrary 
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bodies,  the  pressure  distribution  can  also  be  computed  directly  from 
equations  (l)  and  (2)  (provided  that  f(|)  is  a  continuous  function 
of  x)  by  either  analytical  or  graphical  procedures;  however,  the 
graphical  approach  involves  extensive  computation  and  a  separate  graphical 
integration  for  each  point  considered,  and  the  analytical  approach 
involves  the  approximation  of  the  body  shape  by  a  suitable  expression 
and  (except  in  special  cases)  difficult  and  tedious  integration. 

Equation  (l)  may  be  written  as 


i=n 

=1 


2Qi 


i=r,\Ax  -  *i)2  -  (p®z)2 


(3) 


where  each  of  the  n  point  sources  located  on  the  axis  of  the  body 
ahead  of  (x  —  PRX)  represents  a  section  of  the  body  of  length  Ax. 
Comparison  of  equation  (3)  with  equation  (l)  shows  that 


Q  =  f»(5)Ax 


(4) 


It  is  evident  that,  if  the  arbitrary  body  may  be  adequately  represented 
by  a  relatively  small  number  of  point  sources,  the  pressure  distribution 
may  be  calculated  rapidly  by  use  of  equation  (3).  The  computation  is 
greatly  facilitated  by  the  direct  proportionality  between  the  pressure 
coefficient  and  the  value  of  Q  for  a  single  source.  (See  equation  (3) 
for  n  =  1.)  This  direct  proportionality  allows  the  pressure  field  for 
a  single  source  to  be  computed  once  as  a  function  of  0r  and  used  for 
any  source  or  sink;  at  any  Mach  number  within  the  limitations  of  the 
basic  equations.  This  universal  pressure  field  may  be  conveniently 
plotted  as  a  ,,grid,,  of  lines  of  constant  pressure  coefficient  on 
the  x,3r  plane.  Values  for  this  grid  computed  from  equation  (3) 

(for  n  =  1,  Q  =  0.01,  5  =  0;  that  is,  a  single  source  having  a  convenient 
value  of  Q  and  located  at  the  origin  of  the  grid)  are  presented 
as  table  I.  The  grid  is  plotted  in  figure  1.  The  pressure  distribution 
over  the  body  due  to  a  single  source  representing  any  section  of  the 
body  of  length  Ax  may  be  readily  obtained  graphically  by  superimposing 
the  origin  of  the  grid  at  the  source  location  on  a  drawing  of  the  body 
shape  (plotted  as  0R  against  x),  and  reading  the  pressure  shown  by 
the  grid  at  a  number  of  points  on  the  body  surface,  and  then  multiplying 
these  pressures  by  the  ratio  of  the  Q  for  the  source  being  considered 
to  the  value  of  Q  for  which  the  grid  was  prepared.  The  process  is 
repeated  for  each  source  necessary  to  represent  the  body,  and  the  pressure 
coefficient  at  any  point  is  obtained  by  adding  the  pressures  at  the 
point  due  to  each  source. 

If  the  body  is  to  be  satisfactorily  represented  by  a  small  number 
of  discrete  sources,  the  location  of  a  source  within  its  section  of  the 
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"body  -will  be  Important.  It  was  found  (see  appendix)  that  the  point  source 
should  be  located  on  the  axis  Ax/4  ahead  of  the  rear  end  of  the  section 
of  the  body  of  length  Ax  which  it  represented  and  the  Q  for  the 
source  should  be  computed  from  the  radius  and  slope  of  the  body  surface 
at  the  source  location.  With  this  source  location  and  strength,  the 
pressures  on  the  body  surface  at  and  behind  the  intersection  of  the  Mach 
line  from  the  axis  at  the  end  of  the  interval  with  the  body  surface 
closely  approximated  the  pressures  obtained  from  a  continuous  distribution 
of  sources  along  the  axis  of  the  section.  As  the  pressures  on  the  body 
due  to  a  number  of  sources  are  summed  only  at  the  intersection  of  the 
Mach  line  from  the  end  of  each  interval  with  the  body  surface  (see 
appendix),  the  infinite  pressures  which  equation  (3)  gives  for  x  =  f-  +  (3r 
(along  the  Mach  line  from  each  source)  are  avoided.  The  validity  of 
results  obtained  by  the  point-source  method  and  the  computational 
procedures  used  are  demonstrated  by  an  example. 


EXAMPLE 


Computation  of  a  pressure  distribution  by  the  point— source  method 
is  accomplished  in  the  following  steps: 

(1)  Divide  selected  body  shape  into  sections  and  determine  source 
locations. 

(2)  Graphically  compute  pressures  on  body  due  to  sources  of  Q  =  0.01 
by  successive  superposition  of  grid  on  body  shape  plotted  as  0E  against 

(3)  Prepare  tabular  form  for  recording  values. 

(4)  Determine  value  of  Q  for  each  source. 

(5)  Compute  pressure  on  body  due  to  each  source. 

(6)  Sum  pressures  due  to  each  source  to  obtain  desired  pressure 
distribution. 

A  parabolic  body  of  fineness  ratio  12  was  chosen  for  this  example 
as  representative  of  bodies  suitable  for  treatment  by  the  linearized 
theory  and  because  it  could  be  easily  treated  by  the  method  of  Jones 
and  Margolis  (reference  4)  for  comparison.  The  pressure  distribution 
was  computed  at  a  Mach  number  of  1.4l. 

Body.—  The  body  shape  is  specified  by 

®  =  “r(l  -  ^2)  (5) 

Equation  (5)  describes  a  parabolic  body  of  unit  length  symmetrical  about 
its  maximum  diameter  which  is  located  at  x  =  0.  The  body  was  arbitrarily 
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divided  into  20  sections  of  equal  length  (Ax  =  0.05),  and,  as  described 
in  the  preceding  section,  a  source  location  was  chosen  Ax/4  ahead  of 
the  end  of  each  section. 

Graphical  computations.—  The  successive  superposition  of  the  grid 
origin  at  the  chosen  values  of  £  is  accomplished  by  plotting  the  grid 
and  body  shape  to-  the  same  scale  on  separate  sheets  of  semitransparent 
coordinate  paper  and  matching  the  two  sheets  on  a  glass  surface  lighted 
from  behind.  A  convenient  size  for  the  grid  for  most  purposes  has  been 
found  to  be  10  by  15  inches  with  the  Pr— scale  10  times  the  x— scale. 

The  grid  (fig.  l)  should  be  plotted  carefully  from  the  data  of  table  I 

and  can,  of  course,  be  used  for  all  bodies  and  Mach  numbers.  The  body 
shape  (plotted  in  the  form  PE  against  x)  is  shown  in  figure  2  with 
the  grid  origin  superimposed  on  the  source  |  =  -O.2625  which  represents 
the  section  of  the  body  between  x  =  -0.30  and  x  =  -0.25. 

Table.—  A  convenient  form  for  tabulating  the  values  obtained  in  the 
computation  is  presented  as  table  II.  The  abscissas  of  each  of  the 
20  sources  are  arranged  vertically  in  two  sets  on  the  left  of  the  table 
and  the  values  of  P  used  for  the  grid  are  used  for  column  headings. 

The  values  of  x  at  which  the  lines  of  constant  P  of  the  grid  cross 

the  body  surface  are  recorded  in  the  appropriate  column  in  the  upper  part 
of  the  table.  For  the  case  shown  in  figure  2  (I  =  — O.2625)  the  values 
of  x  in  table  II  are  underlined. 

The  last  column  in  the  upper  part  of  table  II,  headed  x2,  gives  the 
values  of  x  at  which  the  pressures  due  to  each  source  are  to  be  summed. 
These  values  are  obtained  by  superimposing  the  grid  origin  at  the  end 
of  the  interval  and  reading  the  value  of  x  on  the  body  where  the 
line  of  P  =  co  (Mach  line)  crosses  the  body  surface.  For  the  interval 
just  mentioned  (g  =  -O.2625),  the  grid  origin  is  set  at  x  =  -0.25  and 
the  value  found  for  x^  is  -0.216. 

Determination  of  Q  for  each  source.—  In  this  example  the  values 
of  Q  for  each  of  the  20  sources  may  be  found  easily  by  use  of 
equations  (2),  (4),  and  (5).  In  the  more  general  case  of  the  arbitrary 
body  with  which  this  paper  is  concerned,  however,  the  value  of  Q,  for 
each  source  must  be  determined  from  the  ordinates  of  the  body  by  use 
of  the  differential  forms  of  equations  (2)  and  (4)  as  shown  in  figure  3* 
Any  convenient  increment  may  be  used  to  determine  the  variation 
of  B(dE/dx)  with  x  from  the  variation  of  E  with  x  (lower  part 
of  fig.  (3))*  However,  as  the  value  of  Q  must  be  determined  at  the 
source  location,  A[E(dE/dx)]  is  taken  from  the  midpoint  to  the  end  of 
the  interval  and  multiplied  by  2  (upper  part  of  fig.  3)*  The  ratios 
of  the  value  of  Q  determined  for  each  source  to  the  value  of  Q  used 
to  prepare  the  grid  (0.01)  are  tabulated  in  the  last  column  of  the  lower 
part  of  table  II  opposite  the  value  of  £  for  each  source. 

Computation  of  pressures  on  body  due  to  each  source.—  The  data  in  the 
upper  part  of  table  II  is  the  pressure  distribution  over  the  body  if  all 
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the  sources  have  Q  =  0.01.  As  the  pressure  at  any  point  due  to  one 
source  is  directly  proportional  to  the  value  of  Q,  the  actual  pressure 
distribution  on  the  body  due  to  the  source  is  determined  by  multiplying 


Q 


for  the 


is  entered  in  the 
i  =  -0.2625,  the 


x  =  —0.210 


the  value  of  P  heading  each  column  by  the  ratio 

Qgrid 

particular  source  being  considered.  The  resulting  P 
lower  part  of  table  II.  For  example,  for  the  source 

Q* 

ratio  - —  is  -0.0118  and  the  pressure  coefficient  at 

^grid 

(upper  part  of  table,  line  |  =  —  O.2625,  column  P  =  0.5)  is 
P  =  -0.0118(0.5)  =  -0 . 0059  (lower  part  of  table,  line  |  =  -0.2625, 
column  P  =  O.5).*  For  the  same  source  the  pressure  coefficient 
at  x  =  -0.132  (upper  part,  line  |  =  -0.2625,  column  P  =  0.l6) 
is  P  =  -0.0118(0.16)  =  -0.0019  (lower  part,  line  £  =  -O.2625, 
column  P  =s  0.16)  .  These  values  are  bracketed  in  table  II. 


Summation  of  pressures  due  to  each  source .—  The  pressure  coefficients 
for  each  source  tabulated  in  the  lower  part  of  table  II  are  plotted  against 
the  corresponding  values  of  x  at  which  they  occur  (upper  part  of  table  II) 
in  figure  4.  The  logarithmic  ordinate  scale  is  used  so  that  the  points 
are  separated  at  low  values  of  the  pressure  coefficient.  Curves  are 
faired  through  the  points  for  each  source  and  all  of  the  curves  are  added 
at  the  summation  points  previously  determined  (last  column,  upper  part 
of  table  II).  These  sums  and  the  values  of  x^  at  which  they  occur  are 
the  desired  pressure  distribution  on  the  body.  The  pressure  coefficient 
at  the  nose  of  the  body  is  computed  directly  from  equation  (Al) 

(see  appendix)  under  the  assumption  that  the  extreme  nose  of  the  body 
approximates  a  cone. 

Computation  for  different  Mach  numbers.—  In  order  to  compute  the 
pressure  distribution  on  the  same  body  at  a  different  Mach  number,  the 
body  shape  (in  the  form  PE  against  x)  must  be  replotted  and  the  data 
in  the  upper  half  of  table  II  redetermined.  The  data  in  the  lower  half 
of  table  II  are  not  changed  and  the  pressure  distribution  is  computed 
from  the  new  table  exactly  as  in  the  preceding  example. 

Time  required  for  computation.—  Once  the  grid  has  been  prepared, 
computation  of  a  20-point  pressure  distribution  for  an  arbitrary  body 
by  this  method  required  about  8  man-hours  and,  for  the  same  body,  results 
may  be  obtained  for  a  different  Mach  number  in  about  4  additional  man¬ 
hours.  For  parabolic  bodies  only  the  same  number  of  points  can  be 
computed  by  the  method  of  Jones  and  Margolis  (reference  4)  in  somewhat 
less  than  half  this  time. 
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*■ 

EVALUATION  OF  METHOD 


Comparison  with  results  by  method  of  Jones  and  Margolis The 
pressure  distribution  determined  in  the  example  (fig.  4)  is  reproduced 
in  figure  5  together  with  the  pressure  distribution  computed  for  the 
same  body  by  the  method  of  Jones  and  Margolis  (reference  4) .  The 
expression  for  this  pressure  distribution  may  be  obtained  from 
equations  (l),  (2),  and  (5)  and  is  (see  reference  4} 


1 

+  6  [(x  -  xq)  -  4xj \J (x  —  Xq)2  -  (PEX)2  ’ 


(6) 


Examination  of  figure  5  shows  that  the  pressure  distribution  determined 
by  the  present  method  agrees  with  that  of  Jones  and  Margolis  except  for 
a  small  and  apparently  random  scatter  which  approximately  corresponds 
to  the  uncertainty  of  the  graphical  and  numerical  computations. 


In  order  to  investigate  further  the  uncertainties  introduced  by  the 
approximations  involved  in  the  point-nsource  method,  the  pressure  distri¬ 
bution  was  calculated  for  a  parabolic  body  of  fineness  ratio  6. 

Inasmuch  as  large  pressure  coefficients  and  gradients  are  encountered 
for  the  body  of  fineness  ratio  6  at  M  =  1.4l,  discrepancies  should 
be  more  easily  discernable  than  in  the  case  of  the  body  of  fineness 
ratio  12.  The  results  of  this  calculation  for  10  and  20  intervals 
(Ax  =  0.10  and  0.05)  are  compared  in  figure  6  with  results  calculated 
from  equation  (6) .  This  fineness  ratio  was  chosen  because  it  approxi¬ 
mates  the  lower  limit  of  fineness  ratio  for  which  the  linearized  theory 
gives  a  completely  supersonic  flow  field  about  the  body  at  a  Mach 
number  of  1.4l.  It  is  interesting  to  note,  however,  that  recent  results 
obtained  by  the  MCA  wing— flow  method  showed  that  pressure  distributions 
of  the  type  predicted  by  Jones  and  Margolis  occurred  for  a  body  of 
fineness  ratio  6  at  a  much  lower  supersonic  Mach  number  even  though  a 
large  subsonic  region  was  present  at  the  nose. 


Figure  6  shows  that  once  more  the  results  of  the  present  method 
'.agree  closely  with  those  of  Jones  and  Margolis.  The  scatter  of  the  data 
in  figure  6  (body  of  fineness  ratio  6)  appears  to  be  somewhat  greater 
than  the  scatter  of  the  data  in  figure  5  (body  of  fineness  ratio  12). 
Such  a  result  might  be  expected  because  of  the  steeper  gradients 
encountered  on  the  body  of  fineness  ratio  6.  Little  difference  is 
evident  in  figure  6  between  the  results  for  the  two  different  intervals. 


9 


MCA  TN  No .  1768 


9 


As  additional  computations  for  the  intervals  of  0.02  (50  sections)  near 
the  nose  of  the  "body  shoved  no  significant  difference  from  the  results 
presented,  it  is  concluded  that  the  results  obtained  are  not  appreciably 
affected  by  the  interval  size  (within  the  range  investigated)  and  that 
a  parabolic  body  may  be  satisfactorily  represented  by  as  few  as  10  sources. 
The  method  does  not  require  that  a  constant  interval  length  be  used  and 
for  arbitrary  bodies  it  is  recommended  that  shorter  intervals  be  used  in 
regions  where  the  body  cross  section  is  changing  rapidly. 

The  small  discrepancies  between  the  results  computed  by  the  point- 
source  method  and  the  method  of  reference  4  are  negligible  in  comparison 
with  the  uncertainties,  introduced  into  both  methods  by  the  assumptions 
involved  in  the  linear  theoiy  from  which  equations  (l)  and  (2)  are  taken. 

The  linear— theory  results  have  been  shown  by  several  authors  (for 
example,  reference  2)  to  be  in  error  by  as  much  as  20  percent  at  some 
Mach  numbers  in  the  case  of  a  10°  cone. 

Comparison  of  results  by  point— source  method  with  results  by 
characteristics  methods  and  with  wind— tunnel  measurements.—  Some 
information  on  the  usefulness  and  reliability  of  the  point-source  method 
may  be  obtained  by  comparison  of  results  with  results  of  calculations 
by  the  characteristics  method  and  with  wind-tunnel  measurements  for 
the  German  missile  A4VTP  presented  in  references  7  and  8.  Pressure 
distributions  measured  at  Mach  numbers  of  I.56  and  I.87  are  presented 
in  reference  7  and  these  measurements  are  compared  with  calculations  by 
the  graphical  characteristics  method  of  Tollmein  and  the  computational 
characteristics  method  of  Sauer  in  reference  8.  The  body  shape  of 
the  A4V1P  is  not  particularly  suited  for  calculation  by  linearized 
methods  because  of  the  bluntness  of  the  nose  (ogival  nose  having  a  nose- 
length  to  maximum  diameter  ratio  of  3*5*  nose  half  angle  15-9°-  However, 
the  A4V1P  is  the  only  body  for  which  relatively  complete  experimental 
data  are  available,  and  the  comparison  should  indicate  the  accuracy  with 
which  the  method  predicts  the  pressure  distribution  on  such  relatively 
blunt  bodies. 

The  measured  pressure— distribution  points  and  fairing  (solid  line) 
reproduced  from  reference  7  are  shown  in  figure  7  together  with  the 
results  obtained  by  the  method  presented  herein  (dash  line).  Examination 
of  this  figure  shows  that  the  point-source  method  satisfactorily  reproduces 
the  salient  features  of  the  measured  pressure  distribution  but  shifts  the 
sharp  breaks  which  occur  at  the  beginning  and  end  of  the  cylindrical 
section  toward  the  rear  of  the  body.  The  magnitude  of  the  maximum  positive 
and  negative  pressures  are  also  somewhat  ^exaggerated. 

The  fairing  of  the  experimental  data  and  the  pressure  distribution 
calculated  by  the  point-source  method  shown  in  figure  7  nre  reproduced 
in  figure  8  where  they  are  compared  with  calculations  by  the  characteristics 
methods  of  Sauer  and  Tollmein  (reference  8).  Figure  8  shows  that  Sauer1  s 
calculations  agree  closely  with  the  experimental  results  and  Tollmein1  s 
agree  less  closely,  particularly  near  the  middle  of  the  nose  and  on  the 
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cylindrical  section  of  the  body.  The  calculation  by  the  point— source 
method  agrees  closely  with  that  of  Tollmein  near  the  nose  of  the  body,  the 
most  obvious  discrepancy  being  the  exaggeration  and  rearward  shift  of 
the  minimum  pressure  peak  near  the  beginning  of  the  cylindrical  section 
previously  mentioned. 

A  comparison  similar  to  that  shown  in  figure  8  is  presented  in 
figure  9  for  a  Mach  number  of  1.5 6.  In  general,  this  comparison  shows 
the  same  results  as  does  figure  8;  however,  the  overestimation  of  the 
pressures  on  the  nose  by  the  characteristics  methods,  which  was  barely 
visible  in  figure  8,  is  increased  appreciably  at  the  lower  Mach  number. 

It  will  be  noted  that  none  of  the  calculated  pressure  distributions 
agree  with  the  measurements  on  the  rear  ogive  of  the  missile.  Reference  7 
states  that  the  pressures  in  this  region  are  greatly  affected  by  the 
four  large  stabilizing  tail  fins  which  were  in  place  when  the  measurements 
were  made.  Also,  these  pressures  could  be  subject  to  large  boundary- 
layer  effects. 

The  sharp  breaks  in  the  pressure  distributions  for  the  A4V1P  are 
characteristic  of  the  class  of  bodies  having  cylindrical  center  sections 
faired  into  a  curved  nose  and/or  tail  section  so  that  the  second 
derivative  of  the  radius  with  respect  to  x  is  discontinuous.  It  is 
apparent,  however,  that  the  minimum  pressure  peak  computed  from 
equations  (l)  and  (2)  always  occurs  at  a  distance  pR  behind  the 
beginning  of  the  cylindrical  section  rather  than  at  the  beginning  of  the 
cylindrical  section  where  the  experimental  results  and  rigorous  theoretical 
treatment  show  it  must  be  located.  The  occurrence  of  the  rearward  shift 
for  one  class  of  bodies  suggests  that  such  a  shift  may  be  present  in  the 
pressure— distribution  results  computed  from  equations  (l)  and  (2)  for 
other  types  of  bodies.  Extensive  experimental  pressure-distribution 
measurements  are  required  to  verify  the  existence  of  the  shift  for 
arbitrary  bodies,  however. 

Inasmuch  as  the  approximations  involved  in  the  linearized  theory 
become  more  representative  as  the  body  becomes  more  slender  and  because, 
in  general,  slender  bodies  are  better  suited  for  aerodynamic  purposes, 
the  method  presented  herein  Is  believed  to  provide  a  useful  means  for 
rapid  estimation  of  the  pressure  distribution  on  slender,  arbitrary 
bodies  of  revolution  suitable  for  use  as  supersonic  aircraft  or  missiles. 
Insufficient  experimental  evidence  is  available  to  establish  definitely 
the  limits  of  applicability  of  the  method;  however,  it  appears  doubtful 
that  results  of  usable  accuracy  could  be  obtained  for  bodies  more  blunt 
than  the  A^YIP  used  as  an  example  herein. 

Illustration  of  change  in  pressure  distribution  resulting  from  a 
moderate  change  in  body  shape.—  In  order  to  illustrate  the  magnitude  of 
the  changes  in  pressure  distribution  which  result  from  moderate  change 
in  body  shape  the  pressure  distribution  calculated  for  the  test  body  of 
fineness  ratio  12  used  in  the  Langley  Free-Fall  Test  Program  is  shown 
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in  figure  10  compared  with  that  for  a  parabolic  body  of  the  same  fineness 
ratio.  The  body  shapes  are  compared  in  the  lower  part  of  the  figure. 

The  free-fall  test  body  was  derived  from  an  NACA  low-drag  airfoil  and 
had  a  cusped  tail  faired  into  a  tail  boom  and  a  relatively  blunt  nose. 

Some  indication  of  the  over-all  accuracy  of  the  pressure  distribution 
calculated  for  the  free-fall  test  body  may  be  obtained  by  comparison 
of  the  drag  coefficient  obtained  by  integrating  the  pressure  distribution 
with  the  experimentally  measured  value.  The  computed  value  (O.O8O)  is 
in  good  agreement  with  experimental  results  if  a  shirt— friction  coeffi¬ 
cient  of  0.003  based  on  surface  area  is  assumed.  Examination  of 
figure  10  shows  that  the  free-fall  test  body  is  more  critical  than  the 
parabolic  body  with  respect  to  flow  separation  because  of  the  steep 
adverse  gradient  encountered  near  the  rear  of  the  body . 


CONCLUDING  REMARKS 


A  method  is  presented  by  which  the  pressure  distribution  on  an 
arbitrary  body  of  revolution  moving  at  supersonic  speeds  may  be  calculated 
rapidly.  The  method  was  developed  from  the  linearized  potential— flow 
relations  reported  by  Jones  and  Margolis  in  NACA  TN  No.  108l  and  by 
other  authors. 

The  method  is  essentially  a  graphical  computing  technique  based  on 
the  assumption  that  a  body  may  be  represented  satisfactorily  by  a 
relatively  small  number  of  point  sources  and  sinks.  Pressure  distri¬ 
butions  computed  for  parabolic  bodies  by  this  method  were  in  excellent 
agreement  with  the  results  of  Jones  and  Margolis.  As  few  as  10  sources 
m-iri  sinks  satisfactorily  represented  a  parabolic  body  of  fineness 
ratio  6. 

Comparison  of  pressure  distributions  calculated  for  the  German 
missile  AtVlP  with  wind— tunnel .measurements  showed  that  the  point-source 
method  reproduced  the  salient  features  of  the  measured  pressure 
distribution,  but  that  the  minimum  pressure  peak  was  exaggerated  and 
shifted  to  the  rear. 

Pressure  distributions  for  the  AtVIP  at  Mach  numbers  of  1.87  and  I.56 
computed  by  the  point-source  method  were  somewhat  less  satisfactory 
than  those  computed  by  the  characteristics  methods  of  Tollmein  or  Sauer. 
However,  as  the  time  required  to  compute  a  pressure  distribution  by 
the  point-source  method  is  appreciably  less  than  is  required  by  either 
of  the  two  characteristics  methods  and  as  the  salient  features  of  the 
pressure  distribution  are  satisfactorily  reproduced,  the  point-source 
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method  is  believed  to  provide  a  useful  means  for  rapid  estimation  of 
the  pressure  distribution  on  slender,  arbitrary  bodies  of  revolution 
suitable  for  supersonic  aircraft  or  missiles. 
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APPENDIX 


LOCATION  OF  A  SING-LE  SOURCE  EQUIVALENT  TO  A  CONTINUOUS 
DISTRIBUTION  OF  SOURCES  WITHIN  A  FINITE  INTERVAL 


In  the  linearized  theory  of  supersonic  flows  (references  1  to  5) 
todies  of  revolution  are  assumed  to  he  represented  by  a  continuous 
distribution  of  sources  and  sinks  along  the  axis  of  the  body  and  the 
pressures  resulting  from  these  sources  and  sinks  are  assumed  to  be 
propagated  only  within  the  Mach  cones.  Thus,  if  an  interval  of 
length  Ax  extending  from  the  vertex  of  a  semi— inf inite  cone  is 
considered^  it  is  apparent  that  only  at  the  intersection  point  (the 
intersection  of  the  Mach  cone  from  the  .  axis  at  the  end  of  the  interval 
with  the  cone  surface)  is  the  pressure  influenced  by  all  of  and  only  by 
that  part  of  the  continuous  source  distribution  within  the  interval.  As 
the  rate  of  change  of  source  strength  f1  (|)  is  constant  for  a  cone_, 
the  location  of  a  single  source  which  causes  the  same  pressure  at  the 
intersection  point  as  does  the  continuous  distribution  of  sources  within 
the  interval  may  be  determined  in  the  following  manner. 

The  pressure  coefficient  on  the  surface  of  a  cone  of  semivertex 
angle  6  due  to  a  continuous  distribution  of  sources  along  the  axis  may 
be  found  from  equations  (l)  and  (2)  to  be  constant  and  equal  to 

P  =  2tan%  cosh""^^CQ^—- ^  (Al) 

The  pressure  coefficient  on  the  surface  of  the  cone  at  any  point  x 
due  to  a  single  source  located  at  I  is  (see  equation  (3)  for  n  =  l) 


(A2) 


R  =  x  tan  0 


then  from  equations  (2)  and  (4) 


As  for  the  cone  (with  vertex  at  the  origin) 


Q  =  Ax  tan^8 
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Equating  the  pressure  coefficients  given  by  equations  (Al)  and  (A2) 
at  the  intersection  point  and  expressing  the  location  of  the  source  ahead 
of  the  rear  end  of  the  interval  as  a  fraction  of  the  Interval  length 

(&  =  Ax  -  t  =  COp—-^)  gives 


Computations  from  equation  (A3)  for  values  of  t  from  1.01  to  10.0 
show  that  the  ratio  B/Ax  is  approximately  equal  to  0.25.  The  maximum 
variation  in  this  range,  which  corresponds  to  cone  angles  of  5°  to  4-5° 
at  M  =  l,4l  (or  Mach  numbers  Trom  1.15  to  5-75  for  a  10°  cone),  is 
less  than  ±1  percent  of  Ax. 

An  expression  for  the  pressures  on  the  surface  of  the  cone  behind 
the  intersection  point  due  to  the  continuous  distribution  of  sources 
within  the  interval  may  be  obtained  from  equation  (l)  by  use  of  appropriate 
limits  of  integration.  The  variation  of  pressure  with  x  given  by  this 
expression  is  closely  approximated  by  equation  (A2)  for  the  source 
location  just  determined. 

It  is  therefore  evident  that  the  pressure  on  the  surface  of  a  cone 
at  and  behind  the  intersection  point  due  to  a  continuous  distribution 
of  sources  along  the  axis  within  an  interval  Ax  from  the  vertex  is 
closely  approximated  by  a  single  source  located  Ax/4  ahead  of  the  end 
of  the  interval.  Further  computation  along  these  lines  has  indicated 
that  this  principle  may  be  applied  to  an  interval  located  in  any  part 
of  a  conical  or  parabolic  body  of  revolution  with  reasonable  accuracy, 
provided  that  the  pressure  function  Q  is  computed  (equations  (2)  and  (4)) 
from  the  radius  and  slope  of  the  body  surface  at  the  source  location. 

It  appears  reasonable  to  assume,  therefore,  that  within  the  general 
limitations  of  the  linear  theory,  the  principle  may  be  applied  to 
arbitrary  bodies. 
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TABLE  I 

VALUES  FOR  GRID  (FIG.  l)  OF  LINES  OF  CONSTANT  PRESSURE  COEFFICIENT 
IN  x,  Pr  PLANE  CALCULATED  FROM  EQUATION  (3);  SOURCE  AT  ORIGIN; 
PRESSURE  FUNCTION  Q  =  O.Ol 


\Pr 
P \ 

0 

0.01 

0.02 

0.03 

0.05 

0.09 

X 

00 

0 

0.0100 

0.0200 

0.0300 

0.0500 

0.0900 

2.00 

.0100 

.0141 

.0224 

.0316 

.0510 

.0908 

1.00 

.0200 

.0224 

.0283 

.0361 

.0538 

.0922 

.70 

.0286 

.0303 

.0349 

.04l4 

.0576 

.0944 

.50 

.0400 

.0412 

.0447 

.0500 

.0640 

.0995 

.30 

.0667 

.0674 

.O696 

.0731 

.0834 

.1120 

.20 

.1000 

.1005 

.1020 

.1044 

.1118 

.1345 

.16 

.1250 

.1254 

.1266 

.1286 

.1346 

.1540 

.10 

.2000 

.2002 

.2010 

.2022 

.2062 

.2193 

.08 

.2500 

.2502 

.2508 

.2518 

.2550 

.2657 

.05 

.4000 

.4001 

.4005 

.4011 

.4031 

.4100 

•  03 

.6667 

.6667 

.6670 

.6674 

.6685 

.6727 

.02 

1.0000 

1.0000 

1.0002 

1.0004 

1.0012 

1.0040 

MCA  TN  No.  1768 


IT 


TABLE  II 

FORM  FOR  TABULATING-  VALUES  OF  P  ON  SURFACE  OF  BODY  DUE  TO  EACH  SOURCE  AND  VALUES  OF  x  AT  WHICH 
THEY  OCCUR*  EXAMPLE:  PARABOLIC  BODY  OF  FINENESS  RATIO  12  AT  A  MACH  NUMBER  OF  lAl 
[Underlined  and  bracketed  values  are  for  example] 
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Figure  5.-  Comparison  of  pressure  distributions  Figure  6.-  Comparison  of  pressure  distributions 

for  a  parabolic  body  of  fineness  ratio  12  for  a  parabolic  body  of  fineness  ratio  6 

computed  at  a  Mach  number  of  1.41  by  the  computed  at  a  Mach  number  of  1.4l  by  the 

point -source  method  and  by  the  method  of  point -source  method  and  by  the  method  of 

Jones  and  Margolis  (reference  4).  Jones  and  Margolis  (reference  4). 
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Figure  7.-  Comparison  of  calculated  and  experi -  Figure  8.-  Comparison  of  pressure  distributions 

mental  pressure  distributions  for  the  calculated  by  the  present  method  and  by  the 

A4V1P  at  a  Mach  number  of  I.87.  Expert-  characteristics  methods  of  Sauer  and  Tollmein 

mental  points-  and  solid  fairing  taken  with  wind-tunnel  results  for  the  A4V1P  at  a 

from  reference  7.  Mach  number  of  1.87. 
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Figure  9 •  -  Comparison  of  pressure  distributions  Figure  10.-  Comparison  of  pressure  distributions 

calculated  by  the  present  method  and  by  the  calculated  at  a  Mach  number  of  1.4l  by  the 

characteristics  methods  of  Sauer  and  Tollmein  point -source  method  for  a  parabolic  body  of 

■with  wind-tunnel  results  for  the  AkVIP  at  a  fineness  ratio  12  and  for  the  free-fall  test 

Mach  number  of  1.56.  body  of  fineness  ratio  12. 


